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Abstract 

The alternating integrable spin chain and the RSOS(qi, q2',p) model in the presence of a 
quantum impurity are investigated. The boundary free energy due to the impurity is derived, 
the ratios of the corresponding g functions at low and high temperature are specified and their 
relevance to boundary flows in unitary minimal and generalized coset models is discussed. 
Finally, the alternating spin chain with diagonal and non-diagonal integrable boundaries is 
studied, and the corresponding boundary free energy and g functions are derived. 



1 e-mail: doikou@lapp.in2p3.fr 



1 Introduction 



Two dimensional exactly solvable models with boundaries have attracted a great deal of research 
interest recently from the point of view of boundary conformal field theory P , and critical behavior 
[2] , but also because of the rich variety of physical phenomena they display, which in principle can 
be exactly investigated (see e.g. [H] [7|). There has been also much interest on problems related 
to quantum impurities, mainly because of the important role they play in low dimensional physics 
|H], but also because of their relevance to boundary conformal field theory [TJ El El UOj - There are 
numerous studies related to quantum impurities yielding a great number of interesting and useful 
results (see e.g. [TTJ [T3J USD- In this article we focus basically on the thermodynamic analysis of 
lattice integrable systems in the presence of quantum impurities and integrable open boundaries. 

For both integrable lattice models and relativistic integrable field theories, in the bulk, the 
corresponding free energy has been derived and the conformal properties have been extensively 
studied ^1] [21] by means of the thermodynamic Bethe ansatz. It is however of great interest to 
extend these studies for integrable models with boundaries. In analogy to the bulk case, when 
boundaries are added, the corresponding boundary free energy and the so called g function, which 
characterizes the ground state degeneracy due to the boundaries, E3 EE1 EH EH EEj can be 
specified by means of the thermodynamic Bethe ansatz. On the other hand statistical systems 
at the critical point it is known to display conformal invariance |3UJ EH], therefore they can be 
associated with certain conformal field theories. The low temperature behavior of the free energy 
per unit length of such system, in the bulk, is described by [321 EH] 

f(T) = f --T 2 + ..., T«l, (1.1) 
bu 

where c is the central charge of the effective conformal field theory. Furthermore, when boundaries 
are added the obtained free energy is modified up to an contribution (L denotes the size of the 
system) namely, [2| 

f(T) = f (T)-j\ng, (1.2) 

where the ground state degeneracy g is expected to be related with the boundaries of the system. 
One of the main aims, when studying such systems is to employ proper techniques in order to 
specify the central charge and the ground state degeneracy — when boundaries are present. 

As already mentioned from the integrable systems point of view, the central charge and the 
ground state degeneracy g can be identified by employing thermodynamic Bethe ansatz techniques 
(see e.g. [H] [21, HH1 ESI EH1 EH 123) In this stud Y in particular the alternating integrable 
spin chain [33j and the RSOS(qi, q2',p) model [33] in the presence of a quantum impurity ("Kondo 
type" boundaries see e.g. [H1H2]) are investigated via the thermodynamic Bethe ansatz, and the 
corresponding free energy is derived at low and high temperature. The relevance of the results for 
the RSOS(qi, q2',p) model to the boundary flows in minimal [§] and generalized coset models [TU] 
is discussed. Finally the alternating open spin chain with diagonal and non-diagonal integrable 
boundaries is considered and the corresponding boundary free energy and the g functions for the 
left and right boundaries are determined by first principle calculations at low and high temperature. 
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2 Quantum Impurity 



2.1 The alternating spin chain 

Let us first focus on the alternating spin chain in the presence of a quantum impurity. For 
what follows it is necessary to introduce the basic constructing element of the model, namely the 
R matrix, which is a solution of the Yang-Baxter equation HZj 

-Ri2(Ai — A 2 ) i?i3(Ai) .^23(^2) = -^23(^2) -^13(^1) Rn(^i — A 2 ) . (2.1) 
We consider the R matrix obtained in 38. , namely 

R 1 ^ = sinh/i(A + i(~ + a 3 <g> S 3 )) + sinh/i(<7+ ® S~ + a~ <g> S + ) (2.2) 

where S 3 , S act in general, on a 2S + 1 dimensional space V = C 2S+l , and they satisfy the 
following commutation relations 



+ -I _ sinh2z fiS 3 
J smh i\x 



-- ±S ± , 

J smht/i L J 

S 3 \0) = S\0), S + |0> = 0. (2.3) 

We can now define the transfer matrix of the chain 2 

t = tr T (X) (2.4) 

where 

T (A) = i?02JV+l(^ - ®) R 02NW R 02N~lW ■ ■ ■ ^crcM^oiM > ( 2 - 5 ) 

and R l is related to the spin Si = % {% — 1,2) representation (i? 9 is related to the spin | 
representation) (|2.2|) . Following the standard Bethe ansatz method described in e.g. jHlj, |39]-|4T|. 
the Bethe equations are obtained 

M 

e qi (X a ) N e q2 (X a ) N e q (X a - 9) = - JJ e 2 (A Q - A/3) (2.6) 

/3=1 

where 

sinh/i(A + f ) 7T 
6 " (A;;y) = sinhMA-f) ' ~{x (2 ' 7) 

Notice that the main difference between the usual bulk case jMj (without impurities) and (|2.6J) 
is the appearance of the e g term in the left hand side of (|2.6j) due to the presence of the spin | 
impurity. Moreover we should mention that qi, q2 — qi play the role of "flavors" in accordance 
with the picture in [T2], where the spin | (/ flavor) chain is studied with impurity of spin s. 

2 If we derive the transfer matrix of the model with proper inhomogeneities we obtain massless relativistic 
dispersion relations for the particle-like excitations of the model 41 . 
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In the thermodynamic limit, N — > oo, the string hypothesis is valid ^Hl El EH1> namely the 
solutions of equation (|2.6|) can be grouped into strings of length n with the same real part and 
equidistant imaginary parts 

A^) = A£ + ~(n + l-2j), j = l,2,...,n, 
A°' s = (2.8) 

where A™ is real, and X°' s is the negative parity string and A is real. In order to formulate the 
thermodynamic Bethe ansatz we consider all the strings with length n = 1, . . . , v — 1 plus the 
negative parity string. The Bethe ansatz equations (|2.6|) can be written after we apply the string 
hypothesis [T^J HI] 

2 U-l Mm M 

x nq {K - e) n Xn qj (K) N = c-r n n - a?) n g^k - ^ m 

j=l m=l 0=1 f3=l 



2 v-1 M m M 

gM - Q)UgM) N = - n n G im(K - a™) n ^ - a°), (2.10) 

i=l m=l /3=1 /3=1 



where 



. . . , %tx . cosh u(A + „ , 
g „(A;,) = e.(A ± -) = co3h W A _|j , (2.11) 

and X nm , E nm , and G nm are given in the appendix (jA.lJl . 

We introduce the densities of the holes p n and pseudo-particles p n , and once we take the 
logarithm and the derivative of the Bethe ansatz equations (|2.9j) . (|2.1()jl we conclude 

1 1 v ~ x 

Pn(X) = ^( Z n qi (X) + Z nq2 (X)) + —Z nq {\ - 9) - 2J ^nm * /3m (A) - B ln * p (X) 
1 L m=l 

-(p (A) + p (A)) = \{b qi {\) + 6 W (A)) + jb q (X - 9) - £ B lm * p m {X) - a 2 * p (X) (2.12) 

m=l 

where po is the density of the negative parity string, and L = 2N is the length of the spin chain 3 , 
and the quantities Z nm , A nm , B nm , and b n are given in the appendix (jA.5|) . ()A.6|) . (jA.7|) . and ()A.4|) . 
It is also convenient to solve p n (A) in terms of p n (A), therefore we consider the convolution of the 
first of the equations (j2.12)) with the inverse of A nm 

Kiln - $nm - s(u)(8 nm+1 + <5 n m-l) (2.13) 

where s(u) — 2co * h g , s(A) = 2cos 1 h7rA . Having in mind the identities 

A-^*Z mqi (X) = s(X)8 nqi , A£*B 1m {\) = -s{\)6 w - 2 (2.14) 

3 We treat the impurity, which sits at the 2N + 1 site of the chain, separately, therefore L = 2N is the length of 
the bulk part. 
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we obtain the following expressions 



Pn(A) = -s{X){8 ngi + 5 nq2 ) - A nm * Pm{X) - 5 nu - 2 S * p (X) 

A m=l 



1 

L 

p (A) = -p (A) + s*p i/ _ 2 (A) (2.15) 



+ 5 ni/ - 2 s * s * p (A) + T s(X - Q)5 nq 



note the extra j contribution term in the first of the two above equations, which is due to the 
impurity. Let us now derive the free energy of the system f = e — Ts, which is given by 



f = e-Ts = -- E / dX(Z nqi (X) + Z n?2 (A)K(A) - - / d\(b qi (X) + b q2 (X))p (X) 
TJ2 r dX( Pn (X) ln(l + ^) + p n (X) ln(l + ^ 

~,J-00 V Or, A p 



m=l 



„=1™ Pn(X) p n (X)' 

T T dX( Po (X)Hl + ^) + p (X)Hl + ^)). (2.16) 
J-oo v p (A) Po(A) J 

The thermodynamic Bethe ansatz equations are obtained by minimizing the free energy (Sf = 0) 
and by virtue of ()2.12j) . we conclude that they coincide with the ones of the model without 
impurities [12] and they are given by 

Tln(l + ^(A)) = -i(Z ngi (A) + Z rig2 (A))+TEAr l m*ln(l + 77^( 

- T5 ln *ln(l + 7 ?0 - 1 (A)^ 

^l + Vo (A) y 2 ^ v 

- Ta 2 * In (l + T^A)), (2.17) 

where rj n = Alternatively we can write the thermodynamic Bethe ansatz equations by virtue 
of (|2.15J) in the following form, r] n (X) = e -2 ^, 

e n (A) = s{X)*Tln{l+ Vn+1 {X)){l + i ln - 1 {X))-^s{X){8 nqi + 8 ng2 ) 

+ tf w _ 2 a(A)*Tln(l+ 9 7o 1 (A)) 
e„_i(A) = S (A)*Tln(l+^_ 2 (A)) 

e (A) = -s(A)*Tln(l+^_ 2 (A)). (2.18) 

Although the thermodynamic Bethe ansatz equations remain the same as in the bulk case, the free 
energy of the model is modified because of the presence of the quantum impurity. In particular, 
there is a non trivial contribution to the free energy which is given by the following expression, 
after we apply (12~T5T) and (I2~T7I) . (I2~TH1) to (l2~ToT) 

f = e Q + fo + fb + 0(^) (2.19) 



where /q is the bulk free energy given by 

7 1 POO 

fa = -- / dXs(X) ln(l + % (A))(1 + Tfc(A)). (2.20) 

Z J — oo 

eo is the energy of the the state with the seas of strings with length qi, q 2 filled 

eo = -t E / dA^(A) a (A) -rrE d\Z qq .{\ - 9)s(A) (2.21) 

where the non-trivial -| contribution to the bulk ground state energy is due to the impurity. 
Finally, /& is the free energy contribution of the impurity (after we make the shift A — > A — - In T) 

fb = ~ [°° d\s(\--\nT)\n(l + r) q (X)). (2.22) 

We can accurately evaluate differences of boundary free energies, or ratios of g functions (jl.2j) . for 
different temperatures and not specific values at each temperature. This happens basically because 
an overall \ contribution, which can not be explicitly evaluated, may survive in the bulk calculation 
(see ()2.19j0 of the free energy as well (see also [221 1213 121] )• The boundary free energy, and the 
corresponding (^-function ()1.2|) can be computed at any temperature by employing numerical 
methods, however it is possible to make analytical calculations at T — » oo "weak-coupling" and 
T — > "strong-coupling" point (see e.g. (21II21)- 

In particular, for T — > oo the main contribution to the integral in (|2.22j) comes from the A — > oo 
behavior. The corresponding behavior of f] q (X) is given by (see also [H1I12]) 

(1 + Cr 1 = (^T^ n=l,...u-2 
(1 + V^-i)- 1 = K (l + O-^l-i (2.23) 

which obviously does not depend on qi, q 2 . The boundary free energy contribution is for q < v — 1 
"weak-coupling" point (see e.g. 01121), 

/« = _|ln(g + l). (2.24) 
In the special value where q = v — 1 the boundary energy becomes following (|2.23J) 

/r = -^lnM, (2-25) 

which is the half of the expected value. In the isotropic case v oo, q can take any value from 1 
to oo and the impurity contribution is given by (|2.24jh 

When T — > the main contribution to the integral in ()2.22j) comes from the A — > — oo behavior 
- recall that we considered the shift A — > A — - In T. The quantity l + rj n for A — > — oo and T — » 
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is given by E2] 



sin 2/ 



blU ^ gi +2 ^ 

Mir' T 



[l + Vn)' 1 = 7 ~ I TTo ' n = g2 + l,...,i/-2 



(n - g 2 + I) 5 

(l+^-i)" 1 = I, (l + ijgr^l-i (2-26) 

where v = v — q 2 . The situation is more complicated now, because the behavior of ?7 ? (A) depends 
clearly on the flavors (see (|2.26j0 . In particular, for q = the boundary entropy is oc T 2 , and 
this is the completely screened case (see e.g [12J. For q < q 1 

rp sin (*iS+H) 

/0 = -iln y i+2 / (2.27) 



sm 



,91+2, 

this is the behavior of the non-trivial "strong-coupling" point (see also [21 E]) with flavor q 2 and 
impurity spin q. For q\ < q < q 2 

T sin f Zfal±i) ) 

/o = _j_ ln y 2 -; i+2 / , (2.28) 

L sin 



,<j2-<?i+2, 

this case also corresponds to a "strong-coupling" point, with flavor qi — q 2 and a reduced spin 
impurity q — q 2 . For q > q 2 

f b =- 1 -\n(q-q 2 -l) (2.29) 

this is the partially screened case with reduced impurity spin q — q\ . Finally, for q = v — 1 ()2.26|) 

f b =-I-\n{v-q 2 ). (2.30) 

We should emphasize again that the boundary free energy is calculated up to an overall \ contri- 
bution which we are not able to derive explicitly. Therefore, we consider only differences of the 
free energy for different temperatures, and from that via relation (jl.2j) we deduce the ratio of the 
g function, 



,OCA 



ln^ = iln^±% (2.31) 

where ln(l + r)™'°) is derived in (|2.23jh ()2.26|) . The model under study is related to WZW$ q and 
W ZW q2 model, as already concluded in [HI H2j. Indeed the central charge conjectured in |Hj, 
and found in [12] is given by 

- 3qi (2.32) 



qt + 2 5q + 2 
6 



and it is expressed as the sum of the central charges of the 577(2) WZW models at level q 2 and 
Sq (Sq = q2 — qi)- Therefore we expect that our results (|2.24j) - (|2.3U|) . (|2.31|) should be related to 
boundary flows in WZWk models in analogy to the bulk case. Finally, we should note that our 
results for the free energy ()2.24j) - (j2.3Uj) for q\ = q 2 coincide with the ones found in for the spin 
f — % chain with spin s = | impurity. 

In general the impurity sitting at the 2N + 1 site of the chain can be thought as an immobile 
"particle" with constant rapidity 0. Therefore, the particle-like excitations of the chain can 
interact (scatter) with the impurity giving rise to specific scattering amplitudes see e.g. |3E] . 
which we are going to study in detail elsewhere. The scattering should involve in addition to the 
usual XXX part (in the isotropic limit of our model), RSOS type scattering as well [13]. This 
is expected since our chain is related to WZWk model, and it is known (see e.g. [46J) that the 
S'-matrix that describes WZWk models has apart from the SU(2) invariant part an RSOS part 
as well, namely jlH] 

Swzw k = S S u{2) ® Sks OS . (2.33) 

Ssu(2) is the usual XXX S matrix and S^g OS is the S matrix for the RSOS model with restriction 
parameter k + 2 see also 07j . 



2.2 The generalized RSOS(qi, q2',p) model 

It is known that the effective conformal field theory for the critical RSOS(l, 1) model is the unitary 
minimal model A4 U , whereas the critical RSOS(qi,p) model corresponds to the generalized SU(2) 
coset model A4(qi, v — q\ — 2) = su ( 2 )n^^( 2 ^-n- 2 it has been also recently shown jHS] that 
the effective conformal field theory for the generalized critical RSOS(qi,q 2 ;p) model (q 2 > qi) 
consists of two copies of 577(2) coset models, namely M.(qi, v — q% — 2) ® M.[q l) 5q). Our aim 
is to study the boundary behavior of the generalized RSOS(qi,q 2 ;p) model, with "Kondo type" 
boundaries. In particular, the corresponding boundary free energy and the ^-function will be 
derived, and their relevance to boundary flows of conformal field theories CHI will be discussed. 

To describe the model, an orthogonal lattice of 2N + 1 horizontal and M vertical sites is 
considered. The Boltzmann weights associated with every site are defined as 

™(U,-,l m ,yA) = ( \ \ ) ■ (2-34) 

With every face i of the lattice an integer Zj is associated, and every pair of adjacent integers 
satisfies the following restriction conditions |4"5] 

< h+i - h + P < 2P, (a) 

P < +k<2u-P, (6) (2.35) 

where P = qi for i odd, P = q 2 for i even (let q 2 > qi), for % = 1, . . . , 2N and P = q for % = 2N + 1 
for the horizontal pairs, while P = p for the vertical pairs (array type II |1H1)- 
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The fused Boltzmann weights have been derived by Date etal in jSU] and they are given by 

w qi ' 1 (a 1 ,a qi+1 ,b qi+1 ,b 1 \X) ^ ^ JJ w 1 ' 1 (a k , a k+1 , h +1 , b k \X + i(k - ft)) (2.36) 

a2---a qi k=l 

where 62 . . . are arbitrary numbers satisfying |6j — b i+ i\ = 1. ty 1,1 are the Boltzmann weights 
for the SOS(l, 1) model -36|, they are non vanishing as long as the condition (I2.35f a)) is satisfied, 
and for P — 1 they are given by the following expressions 

u;(Z,Z±l,Z,Z=Fl|A) = h(i - A) 



w(Z±1,Z,Z=fU|A) = -Zi(A)- 



Zi 



w(l± 1,1,1 ±l,l\\) = h(wi±X)-± (2.37) 

where, 

h(\) = pQ(\)H(X) (2.38) 
H(X) and 0(A) are Jacobi theta functions and, 

hi = h(wi), wi — Wq + U. (2.39) 

We are interested in the critical case where h(X) becomes a simple hyperbolic function i.e., 

, . , . sinh uA 

h A = — ^, 2.40 
sin /x 

100, p and \x are arbitrary constants. Furthermore, 

p— 2 (ji-i 



>i, 61, 65+1, s+i) = n n - j) + a ) 



fc=0 j=0 

£ f[w qi ' 1 (a k ,b k ,b k+1 ,a k+1 \X + i(k-l)), (2.41) 

0,2— a q k=l 

again 6 2 • • • are arbitrary numbers satisfying |6j — 6 i+1 | = 1, and the pairs ai, a g+ i and 61, 
satisfy (|2.35|) . for P = q. The fused weights satisfy the Yang-Baxter equation in the following 
form 

w pq (a, b, g, f\X)w ps (f, g, d, e\X + fi)w qs (g, b, c, d\fi) 

9 

= ™ qs U\ a, g, e|/iK>, b, c,g\X + pK%, c, d, e| A). (2.42) 
Here we only need the explicit expressions for u^' 1 which are 

•*•'//... 1 /'... 1 /' n\\) . \ ; , - A) 

hi 

,/,t/7 ■■ v , ? ' ? , x , ,.-,.-1, s.u KX + ia) 



(1 + 1,1' + 1,1', l\X)) 


= h qv 


:i( 


-A)/i Q 


(l + l,l'-l,l',l\X)) 


- h Ql ' 


:{( 


-A)/i 6 - 


(l-l,l' + l,l',l\X)) 


= h qi ' 
lb qi- 


:J( 


-A)/i c - 


(1-1,1'- 1,1', l\X)) 


- h qi ~ 

n qi- 


:1( 


-A) /id 



to 

//<■>//_ 1 /'.. 1 /' /I An = /,'/■-! i-\)h h ^ d ~ X ^ 

hi 

l ' / — 1 /' /l.\n -. /, , M /r " A - 

Zi, 



(2.43) 



where 



l + l'-qi , I'-l + qi l-l' + qi , I + 1' + q { ... 

n ' b= n ' C= n ' d= n ' ( 2M ) 



and 

K(X) = lih(X + z(k-j)). (2.45) 

j=0 

It is obvious that w qi,1 (a, b, c, d\X) are periodic functions, because they involve only simple hyper- 
bolic functions fl233), d230|) h(A + iv) = -h(X), v=f), i.e. 

w qi '\a, b, c, d|A + iv) = (-)*w 9i - 1 (a, 6, c, d|A). (2.46) 

Now we can define the transfer matrix of the RSOS(qi, q2',p) model 

2JV-1 

= II «;^ > ^ +1 ,6 i+1 ,6 i |A)ti;^(a i+1> o, +a ,6 i+a ,6 i+1 |A) 

w q > p (a 2N+1 , a 2N+2 , b 2N +2, b 2 N+i\X) (2.47) 

where we impose periodic boundary conditions, i.e. a 2 ^ +2 = a\ and b 2 N+2 — b\. By finding the 
eigenvalues of the transfer matrix we end up with the Bethe ansatz equations of the model, (see 

also ma img) 

M 

LJ-\(X a - e)e qi (X a ) N e g2 (X a ) N = - J] e 2 (X a - \ p ) (2.48) 

/?=i 

It is important to mention that the eigenstates of the model are states which satisfy (see also 

[13113 ESI) 

M= l -{ qi + q 2 )L + ^. (2.49) 

To formulate the thermodynamic Bethe ansatz for the RSOS(q\,q 2 ' 1 p) model we consider all the 
strings ()2.8j) with length n — 1, . . . , v — 2, [13 EE]- The densities of the corresponding holes and 
pseudo-particles are derived from the Bethe ansatz equations f)2.48j) . after we insert all the allowed 
strings, and they satisfy 

p n (X) = hzM(X) + ZM (A)) + jZ$(X - 6) - £ * p m (X)., (2.50) 

1 L m=l 

where L = 2N. From the constraint ()2.49|) it follows that p v - 2 = 0, then the density of the v — 2 
string can be written in terms of the remaining densities 

p,_ 2 (A) = p°(X) - £ afc£m * Pm(X) (2.51) 

m=l 

where a^~ 2 is given in the appendix ()A.4|) with z/ — > i/ — 2, and 

.0/ \ sinhfaf) +sinh(q 2 f) 1 sinh(gf) 

P lWj 4cosh(|)sinh((i/-2)f) L2cosh(|)sinh((i/-2)|)' 1 j 
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By means of the relation (|2.51|) the equation (J2.50|) can be rewritten in the following form 

Pn(A) = \{Z^\X) + Z^\\)) + \z^\X - 0) - £ Afc 3 > * Pm (A), (2.53) 

Z^ 2 \ A^~ 2 ^ are given by ()A.5|) . ()A.6|) . with v — > v — 2. Following the standard procedure 
of minimizing the free energy of the system (Sf = 0, / is derived by (j2.16j0 we obtain the 
thermodynamic Bethe ansatz equations which are the same as in the bulk |35 a 

v-3 

r Z ( "- 2, ('A) + Z l "- 2) (\Y) \- 



Tln (l + r/ n (A)) = -\{Z^\\) + Z^\X)) + £ At n 2) * Tin (l + ^(A)), (2.54) 



m=l 



alternatively we can write 



e n (A) = s(X)*T\n(l + Vn+1 (X))(l + Vn ^(X))-^s(X)(5 rm +5 nq2 ). (2.55) 



The corresponding free energy is 

f(T) = e + f + f b + O(j) (2.56) 

where eo is given by (j2.21|) and fo is the bulk part of the free energy given by ()2.20|) . We wish to 
compute the non-trivial boundary part of the free energy, which is given as in the previous case 
— for the alternating spin chain — by (J2.22j) . Again we can evaluate differences of free energies at 
T — > and T — > oo. For this purpose we need the solutions of 1 + rj q 1)2. 5 5 j) . for T — > oo and 
T — > 0. For T — > oo, the main contribution in ()2.22|) comes for A — > oo [HU EH] the solution is 

sin 2 !'-) 

(1 + Cr 1 - . J^, , n = l,...,i/-3. (2.57) 



sm 



and for T — > the main contribution to the free energy (J2.22|) comes for A — > — oo, and the 
corresponding solution is [35] 



sin 2 ' * 



= . 2 :x + i y - = i,.. (i+o = i 



sin 2 ' 



+ = . 2 Xr + i y n = (ft + l,..., ga -l J + i 

sm — :i V+ 

2/ 



sm 
sm 



+ = TwggW » = ft + l,---^-3. (2-58) 



u-q2 



We observe that the solution at T — * oo does not depend on g^, while at T — > the solution clearly 
depends on the values of Having in mind the above solutions we are ready to derive ratios of 
g functions, in particular 



= I i n ( 1 + V: 



oo \ 
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where 1 + 77^°'° are given by (|2.57|) . ()2.58|) . 

It is worth making some remarks concerning ()2.57|) . ()2.58)1 . and ()2.59|) . In the special case 
where q\ = q 2 = 1 we recover the results of jU] for boundary flows in minimal models A4 u=m+ i, 
(a = q — 1), whereas for q\ = q2 > 1 we recover the results of ^Hj for boundary flows in generalized 
SU{2) coset models M.{qi, v — q\ — 2) = M.{k. I). The corresponding S matrix for the generalized 
coset M.{k,l) model conjectured in jlHI is given by 

S = Srsos ® Srsos> (2.60) 

and in the limit I — ► 00 it reduces to Swzw k given in ()2.33|) . Again from ()2.57|) . ()2.58|) and ()2.59|) 
the structure M.{q\, v — q\ — 2) ® .M(gi,c)g) of the effective conformal field theory is manifest 
(compare with the boundary flow in SU{2) coset models [TUj). 



3 Open Boundaries 

The main aim of this section is the investigation of the thermodynamics of the alternating spin 
chain in the presence of integrable boundaries. To construct the spin chain with boundaries in 
addition to the R matrix another constructing element, the K matrix, is needed. The K matrix 
is a solution of the reflection (boundary Yang-Baxter) equation [3] , 

#i 2 (Ai - A 2 ) tfi(Ai) R 2 i(Xi + A 2 ) K 2 (X 2 ) 
= K 2 (X 2 ) tf 12 (A! + A 2 ) K^) R 21 (X 1 - A 2 ) . (3.1) 

In what follows we are going to use Sklyanin's formalism [3] in order to construct the model with 
boundaries. The corresponding transfer matrix t(X) for the open alternating chain of iV sites and 
Si = S 2 = ^ spins is (see also e.g., |H 34J), 

t(X) = tr iC(A) T (A) # -(A) f (A) , (3.2) 

where 4 

T (A) = Rl N (X)Rl N _ 1 (X) ■ ■ ■ i?o 2 (A)i?Q 1 (A), T (A) = Rl (X)Rl (X) ■ ■ ■ R%_ 10 (X)R m (X), (3.3) 

and K + (X, k + ) = K~(—X — k~)* where ^ , k ± are arbitrary boundary parameters for the 
left and right boundaries, and K~ is the matrix jSJ [7] 

= /sinh^(A + zO Ksinh/i2A \ 
1 ' \ Ksinh/i2A sinh//(-A + iO / ' 1 ' 1 

It is interesting to point out that the diagonal boundaries for the critical XXZ spin chain 
(^i = S 2 — 2) correspond to Dirichlet boundary conditions for the sine-Gordon model |2*5j . 
Presumably the purely anti-diagonal X matrix — which completely breaks the U(l) symmetry — 
should correspond to Neumann boundary conditions for the sine-Gordon model. 

4 We could have considered iV=odd and put an impurity at the N site of the chain. Then we would have an 
impurity contribution to the free energy similar to the one in paragraph 2. 
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3.1 The diagonal K matrix 



We consider first the case in which the K matrix is diagonal, namely k = 0. The corresponding 
Bethe ansatz for the model, are known 01] and they are given by 

M 

e x +{\ a )~ 1 e x -(X a )ei(\ a )g 1 (\ a )e qi (\ a ) N e q2 (\ a ) N \\ e 2 (X a - A /3 )e 2 (A Q + A^) (3.5) 

13=1 

where = 2£ ± ± 1 are boundary parameters for the left right boundaries of the chain and they are 
related to some external magnetic field acting on the boundaries. Here, for simplicity we consider 
x ± to be integers less or equal to u—1. The corresponding densities of pseudo-particles and holes, 
in analogy to ()2.12j) satisfy (again we consider as in the bulk all the strings n = 1, . . . , v — 1 and 
the negative parity string (|2.8|0 

1 »-\ 

p n (X) = Z nqi (X) + Z rm (X) + -K n - A nm * Pm{X) - B ln * p (X) 

L m=l 

1 

- (po(A) + Po(A)) = b qi + b q2 + -K - B im * p m {X) - a 2 * p {X) (3.6) 

L m=l 

where L = N is the length of the chain and, 

K n (u) = a n (u) + b n (u) - Z nx +(u) + Z nx -(u) - 1 

K (tu) = ai(u;) + Si(w)-S a +(w) + 6 a .-(w) + l (3.7) 

Z nm , K, and a n are given in the appendix by (|A.5J) . and ()A.4j) . The unit that appears in the 
expressions for K n , K is a result of the subtraction of a 5(X) term from the densities (see also 
I27j). This subtraction seems necessary for the accurate derivation of the density, because in 
the boundary case A's can take values from to oo, as opposed to the bulk case where A's take all 
the values from — oo to oo. 

The free energy for the model with open boundaries is given by (12.16)1 up to a | factor in front 
of the expression. The appearance of the factor | in front of all the integrals comes from the fact 
that we originally derived the integrals from zero to infinity. After we minimize the free energy, 
we obtain the same thermodynamic Bethe ansatz equations as in (12. 17)1 . ()2.18|) . Finally, by virtue 
of ()2.17|1 . ()2.18j) the following expression for the free energy is obtained 

f = fo + h + 0(j) (3.8) 

where fo is the bulk free energy given by 

T r°° 

f = eo-- dXs(X)\n(l + r ]qi (X))(l + r ]q2 (X)) (3.9) 

/ J — oo 

e is the bulk part of the energy of the state with the seas of strings q±, q 2 filled, 

e o = -T E / dXZ m .(X)s(X). (3.10) 
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and the boundary contribution to the free energy is 

h = "of E / rfA ^( A ) ln (! + ^(A)) + / dXK (\) ln(l + ^(A)). (3.11) 

ZL n _^J — oo J —oo 

From the thermodynamic Bethe ansatz equations (|2.17|) we consider the convolution of s(A) with 
the equation of 1)2.17)1 for n — 1 and n = x^ 1 (with a (-) sign in front of the equation for n = x + ), 
we also consider the convolution of s(A) with the equation for the negative parity string, recall 
also relations 1)2.14)1 . Moreover, we consider the equations 1)2.18)1 for all n = 1, . . . , v — 1 and the 
negative parity string, with a (-) sign infront. We add all the above equations and we end up with 
the following expression for the free energy (at T — > 0, oo) 

T T T 

f b = e b + Sf b = e b + — ln(l + rj v -x) + —s * ln(l + r] x +) - —s * ln(l + rj x -) (3.12) 

where e b is the boundary energy contribution 

e b = [°° d\s(X)(K qi (X) + K l 

'±Ll J —oo ^ 
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1 c oo 1 ^ 

= j /_ oo ^ S (A)(-EK(A) + 6JA) " Z qiX+ + Z qiX -) - 5(X)). (3.13) 
The boundary contribution to the free energy for each boundary then is 

<5/± = ^ ln(l + 77._i) ± ^ Ml + ffe±), (3.14) 

and the corresponding ratios of the g + , g~ (g = g + g~) functions for the left and right boundaries 
QOD are 

In 4- 1 l n ( 1 + ^-i) T 1 l n ( 1 + ^) (3 15) 

%o ±_ 2 ln (l + ^ 1 ) T 4 ln (l + r ? S ± )- (3 - 15) 

The last term of the above equation, for x^ integer less than v — 1, corresponds to the quantum 
impurity contribution (up to a — | factor), which has been already computed explicitly for T — > 
oo, 0, ()2.23)) . ()2.26)) . We focus on the first term of ()3.15)) . which we compute again for T — ► oo, 0. 
It is easy to conclude from the solution for T — > oo, ()2.23j) . ()2.26j) 



1 ln (1 + _ 1 

2 (l + 77°_0 2 "I/- ?2 



h ^^^ -ihJL, (3.16) 



The result (13.16)1 is also valid for the spin S = Si = S2 chain. In the isotropic limit the above 
ratio ()3.16j) becomes unit, i.e. there is no difference in the boundary free energy at low and high 
temperature. In the special case where S% = S2 = \ we recover the result found in |2| for the 
sine-Gordon model with boundaries 5 , and also the result of for the XX Z spin chain with 
open boundaries at zero external magnetic field. Note that for x ± — > 00 there is no boundary 
parameter dependence. 

5 In the repulsive regime, i.e. when no bound states "breathers" exist, we can make the following identification 
^ = A + 1 or P 2 = 8(tt — /u) HOE). 
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3.2 The non— diagonal K matrix 



We consider, for the first time in a spin chain model, the more general case of the open alternating 
spin chain with non-diagonal boundaries. The corresponding Bethe ansatz for the XX Z model 
at roots of unity have been recently derived in j3T] . The novelty of the method described in jM] is 
basically that it does not rely on the existence of a reference state "pseudo-vacuum" . It was shown 
then in |51j . that in the case where [i = (v — p + 1) the problem of finding the eigenvalues 
of the transfer matrix (J3.2)) reduces to a set of functional relations which can be written in the 
following compact form, (see also jSHSSHHI]) 



detM[A qi ' q2 (\)} = 



(3.17) 



where 



M[A qi ' q2 (X)] 



A quq2 


_jiim 








~f 91,92 


Kf m 


—Jo 








f 91,92 


\ 11,12 


_j9l,92 














f<21,<?2 

Jp~l 
















f 91,92 
Jp-1 











92 



A 0.1)9- 

_ f 91,92 
Jv 



_ f 9i,92 \ 






f 91,92 
Jp-2 



(3.18) 



where 



and 



A9U92 j 

f 91,92 W = ~ sinh* //(A + iS 1 + U sinh* //(A + iS 2 + i) 

2 2 smh fi(2X + i) 

(sinh fi(X + i£) sinh /i(A — i£) + k 2 sinh 2 2/iA^ 

f9l,92 {X) = f 9l,92 { _ X _ 2i)> f 9l,12 {X) = /(A + ik) 



(3.19) 



(3.20) 



Let now (Qo(A), . . . , Q P (X)) be the null vector of the matrix (|3.18j) with Qk(X) = Q(X + ik) and 



M 



Q(X) = Y[ sinh /i(A — Xj) sinh //(A + Xj 
then the eigenvalues are given by the following expression 



(3.21) 



A**»(A) = /f*"(-A-t) 



Q{X + t) 



+ /o 91 ' 92 (A) 



Q(X 



Q(X) JU " ' Q(X) 

Finally, from the analyticity of the eigenvalues we obtain the Bethe ansatz equations 

sinh fi(X a + |(2£ - 1)) sinh^(A Q - f (2£ + 1)) + k 2 sinh 2 /i(2A Q - f) 
sinh /i(A a - §(2£ - 1)) sinh fi(X a + §(2£ + 1)) + « 2 sinh 2 n{2\ a + i) 



(3.22) 



M 



5 , i(A Q )ei(A Q )e gi (A a ) e,j 2 (A 



II e2 (^« ~ A ( g)e 2 (A Q + A/j). 

(3=1 



(3.23) 
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Note that the boundary parameters at the left and right boundaries have been considered to be 
the same. In fact, the boundary parameters have to be tuned properly (i.e. £~ = £ + and k~ = k + ) 
so that the method described in |5T| can be applied. Note also that in this case the Bethe ansatz 
states have to satisfy the following constraint, 

M = (S l + S 2 )^- 1 -. (3.24) 

The above constraint follows from the asymptotic behavior of the transfer matrix for A — > oo, in 
particular from ()3.2j) it can be deduced 

t{\ -> oo) ~ _ 2e (^ N +^+ 2 ^+^ N (3.25) 

by comparing the later equation with the eigenvalues ()3.22|) we obtain the constraint ([3.24)1 . 

Let us first consider the case with purely anti-diagonal K matrices, namely k becomes very 
big. Then it is obvious that the above Bethe ansatz equations ()3.23|) take the form 

M 

gi{Ky 2 ei{KT 2 gi{K)ei{\ a )e qi (\ a ) N e q2 (\ a ) N = - n e2 (^« ~ A ( g)e 2 (A a + (3.26) 

p=i 

Relation 1)3.24)1 does not impose in this case any extra restriction on the densities as opposed to the 
case of the RSOS model for which a similar relation 1)2.49)1 holds true, and it imposes restrictions 
on the hole densities (i.e. p v -i = 0). Therefore, the thermodynamic Bethe ansatz equations are 
the same as in the bulk case ()2.17|) . and by following exactly the same procedure as in the diagonal 
case, but now with 

K n {u) = —b n (u) - a n (u) - 1, K (u) = -b\{u) - ax{u) + 1 (3.27) 
we find that the non-trivial boundary part of the free energy becomes (at T — > 0, oo) 

*tf= = ^ln(l + f&). (3.28) 
Moreover the boundary contribution to the ground state energy is given by 

1 c OO 1 ^ 

e b = -^l / ^(^(-^(-MA) - b gi (X)) - 5(A)). (3.29) 

It is worth noticing that the state with the seas of strings with length qi, q2 filled, is an allowed 
state for the case with anti-diagonal boundaries, because the constraint (J3.24)) is satisfied. The 
ratios of the g + , g~ functions for the left and right boundary ()1.2)1 are given by 



gt 4 11 (1 + 4) 



where the function (1 + ^i' 00 ) -1 is given, for T — ► oo and T — > 0, by ()2.23j) . ()2.26jl . Notice that 
for qi = 1 (at T — > 0) 5f^ oc T 2 and this is — as in the presence of an impurity at the end of the 
chain — the completely screened case. 
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We are not going to treat the general case in full detail here, since this is mostly a problem of 
mathematical manipulations, he hope though to report on this in detail elsewhere. Let us however 
present the general concept of such computation for the special simple case where £ = 0, the so 
called "free boundary conditions" jH], then the Bethe ansatz equations become 

M 

e 2 c+i(^a)~ 1 e 2 c-i(A a )e ] ; 2 (A Q )ei(A a )5fi(A Q )e (2l (A a ) Ar e g2 (A Q ) Ar J} ^(^q ~ -V)e 2 (A Q + A^).(3.31) 

/3=i 

It is evident that the later Bethe ansatz equations have a similar structure with ()3.5|) — up to the 
ej~ 2 term in the LHS of ()3.31|) — with £ + = £~ = ( and cosh 2 = —j^ 6 . The boundary part of 
the free energy is given again by (j3.11|) and the quantities K n and K$ are given by, 

K n (u) = -a n (u) +b n (u) - Z nx +{u) + Z nx -(u) - 1 

K (u) = a x (u) -b^u) -b x +(uj) + b x -(u) + 1 (3.32) 

with x ± = 2( ± 1. It remains to treat expression ()3.11|) at high and low temperature, but this 
is relatively easy to do because the solutions for r] n , r] are known at T — > oo and T — > ()2.23|) . 
(|2.26p . Basically we have to focus on the part of expression (|3.11|) that involves the terms ±Z nx ^, 
ifej-T for x non-integer. The remaining x independent part of the boundary free energy (J3.11|) . 
denoted as / 6 , can be easily computed by repeating the steps of the previous section, and it is 
given by (at T — > 0, oo) 

1 r°° 1 ^ T 1 T 1 

= ~2L J-oo dAs(A) (2 SH^W + M*)) - *( A )) + I ln ( X + ^-0 + ln(1 + Vl) - (3 ' 33) 

The x ± dependent part of the boundary free energy, denoted as / 6 (with terms that involve 
±Z n3 .=F, ±6^), which needs to be computed explicitly for x ± non-integers, has the form 

ti 2) = -«fE dX(Z nx -(X)-Z nx+ (\))\n(l + V -\X)) 

Alj n=l J -°° 

+ 2L j /_ oo ^-(A)-6, + (A))ln(l + r ?0 - 1 (A)). (3.34) 

Again once we make the shift A — > A — - InT we can compute / 6 analytically, at T — > 0, oo having 
in mind the expressions for Z nm , b n f|A.2|) . (|A.3p . and 7/ n , 7y at high and low temperatures ()2.23|) . 

______ ————— (_) ^2) f2) 

(|2.2f)jl . More specifically, expression ()3.34j) can be written as / b = / 6 (a; - ) + / 6 (x + ) where at 
T — > 0, cx), 

/fV) = Z nK± (0)ln(l + r ? - 1 ) =F ^ ± ( )ln(l +r ?0 - 1 ). (3.35) 

Naturally, the case is similar when diagonal boundaries with x ± = 2£ ± 1 non-integer, are consid- 
ered, i.e. the x ± dependent part of the free energy ()3.11|) is given again by (j3.35j) . Eventually the 

6 Following |H] we have to chose for the "free boundary conditions" Q = which is not an integer 
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computation of the expression (|3.35|) reduces to a simpler problem, which is the derivation of the 
exact Fourier transforms of Z nx and b x when x is non-integer. In particular, expression (|A.4|) for 
b x is valid as long as x < u, while expression (|A.5|) for Z nx is also valid as long as x > n. 

Let us point out that the results of this section are novel not only for the fused model under 
study, but also for the case q\ = q 2 = 1, which corresponds to the XX Z model (sine-Gordon 
model). 

4 Discussion 

In this investigation we focused on the "boundary properties" of the alternating spin chain and the 
RSOS(q\, q 2 ;p) model, therefore we considered both models with certain integrable "boundaries". 
A quantum impurity was added at the last site of the chain ("Kondo type" boundary see also 
[TT] H2j). and the immediate result was a non-trivial contribution to the ground state energy 
of the system, as well as in the free energy. We were able to explicitly evaluate the non-trivial 
contribution to the boundary free energy at low and high temperature. A similar investigation was 
realized for the RSOS(qi, q 2 ;p) model (see also [S3]) and the obtained results, for special values 
of qi, q2, compared with the known boundary flows in unitary minimal models jH] and generalized 
SU{2) coset theories [TOJ . For general values of qi, q 2 our results are rather novel for both the 
alternating spin chain and the RSOS(qi,q 2 ;p) model. 

Furthermore, the alternating chain with diagonal and non-diagonal boundaries was investi- 
gated, and again the presence of the boundaries resulted in a non-trivial contribution to the free 
energy. This contribution gave rise to the g function ( "ground state degeneracy" ) along the lines 
described in . We were able to derive ratios of the g function for the left and right boundaries, 
at low and high temperatures. In the diagonal case for S% — S 2 — | the boundary parameter 
independent part of the ratio of the g functions coincides with the one found in j2E] and j2H] for 
zero external magnetic field. We believe the results of this section are also novel, in particular for 
the anti-diagonal case our results are rather new even for qi — q 2 — 1. 

For completeness the thermodynamics of the RSOS(qi, q 2 ;p) model with open boundaries [33] 
|33J should be also investigated. It would be also of great interest to extend the thermodynamic 
analysis for spin chains and RSOS models related to higher rank algebras. Furthermore, the 
derivation of the Bethe ansatz equations for non-diagonal boundaries |3T], is an essential step 
towards the investigation of exact non-diagonal boundary S matrices in the spin chain framework. 
Let us finally note that a new non-trivial "dynamical" solution of the reflection equation has been 
recently found [3JI], and it has been realized in the context of the sine-Gordon model. It is a 
challenging problem to apply this "dynamical" solution to the XX Z spin chain at roots of unity, 
and investigate the thermodynamics and the scattering process in this framework. We hope to 
address these questions soon in a future work [37] . 
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A Appendix 

We give the explicit expressions of the functions that appear in the Bethe asnatz equations once 
we apply the string hypothesis, namely 

A nm (A) 6| n _ m -(-i| (A) e , 

n— m+3| 

(A)... 

C(n+m— 3) 

(A)e 

(n+m— 1) (A) 

E nm (\) = e| n _ m |(A)e^_ m+2 |(A) . . .e^ +m _ 2) (A)e (n+m )(A) 

G nm (X) = g\ n -m\(X)gf n _ m+2 \(x) . . .gf n+m _ 2 )(X)g( n+m )(\). (A.l) 

Moreover, 

fln(A) = -^-^-lne n (A), 6 n (A) = \ng n (X), (A.2) 

Z7T aX ztc aX 

i d 

(Z nm (x), A nm (X), B nm {x)) = - — — ln(X„ m (A), E nm (X), G nm (x)). (A. 3) 

Z7T aX 

We finally give the following useful Fourier transforms 

sinh((z/-n)f) A sinh (?f) 

= . \ 7— x n<2u, b n (u) = - — } { , n<v, (A.4) 

smh (if) smh 



sinhf (z/ — max(n, m))^j sinhf (min(n, m))- 2 ^ 
sinh^J smh(~J 



2 coth(^) sinhf (v — maxfn, m))% ) sinhf (minfn, m))^ ) 

2 - — ^ sinh(¥) 2j 1 d 

2 cothm sinhf r^.) s i n hf ^) 

References 

[1] J.L. Cardy, Nucl. Phys. B234 (1989) 581. 

[2] I. Affleck and A.W.W. Ludvig, Phys. Rev. Lett. 67 (1991) 161; I. Affleck, M. Oshikawa and 
H. Saleur, [cowJ-mott/ '980 ^11% (1998). 



18 



[3] I.V. Cherednik, Theor. Math. Phys. 61 (1984) 977. 

[4] E.K. Sklyanin, J. Phys. A21 (1988) 2375; P.P. Kulish and E.K. Sklyanin, J. Phys. A24 (1991) 
L435. 

[5] A. Fring and L. Kdberle, Nucl. Phys. B421 (1994) 159; Nucl. Phys. B419 (1994) 647. 

[6] S. Ghoshal and A. B. Zamolodchikov, Int. J. Mod. Phys. A9 (1994) 3841; A9 (1994) 4353. 

[7] H.J. de Vega and A. Gonzalez- Ruiz, J. Phys. A26 (1993) L519 

[8] C.L. Kane and M.P.A. Fisher, Phys. Rev. Lett. 68 (1992) 1220; E. Sorensen, S. Eggert and 
I. Affleck, J. Phys. A26 (1993) 6757. 

[9] F. Lesage, H. Saleur and P. Simonetti, Phys.Lett. B427 (1998) 85. 

[10] C. Ahn and C. Rim, J. Phys. A32 (1999) 2509. 

[11] V.M. Filyov, A.M. Tsvelik amd PB. Wiegmann, Phys. Lett. 81 A (1981) 175; A.M. Tsvelick 
and PB. Wiegmann, Adv. in Phys. 32 (1983) 453. 

[12] N. Andrei and C. Destri, Phys. Rev. Lett. 52 (1984) 364. 

[13] Y. Wang, Phys. Rev. B60 (1999) 9236; Y. Wang and P. Schlottmann, Phys. Rev. B62 (2000) 
3845; A. P. Tonel, A. Foerster, X.-W. Guan and J. Links, \cond-mat/0 ] 1121l3j 

[14] C.N. Yang and CP Yang, Phys. Rev. 150 (1966) 327; J. Math. Phys. 10 (1969) 1115. 

[15] CP Yang, Phys. Rev. A2 (1970) 154. 

[16] M. Gaudin, Phys. Rev. Lett. 26 (1971) 1301. 

[17] M. Takahashi, Prog. Theor. Phys. 46 (1971) 401; M. Takahashi and M. Suzuki, Prog. Theor. 
Phys. 48 (1972) 2187; M. Takahashi, Thermodynamics of One-Dimensional Solvable Models 
(Cambridge University Press, 1999). 

[18] J.D. Johnson and B.M. McCoy, Phys. Rev. A6 (1972) 1613. 

[19] H. Babujian, Nucl. Phys. B215 (1983) 317; H. Babujian and A. Tsvelik, Nucl. Phys. B265 
(1986) 24. 

[20] L. Mezincescu and R.I. Nepomechie, UMTG-170 (1992). 

[21] Al.B. Zamolodchikov, Nucl. Phys. B342 (1990) 695; Al.B. Zamolodchikov, Phys. Lett. B253 
(1991) 391. 

[22] A.B. Zamolodchikov, Nucl. Phys. B358 (1991) 497, 524; Nucl. Phys. B366 (1991) 122. 
[23] T.R. Klassen and E. Melzer, Nucl. Phys. B338 (1990) 485. 



19 



[24] C. Destri, H.J. de Vega, Nucl. Phys. B438 (1995). 

[25] A. LeClair, G. Mussardo, H. Saleur and S. Skorik , Nucl. Phys. B453 (1995) 581. 

[26] P. Fendley, H. Saleur and N.P Warner, Nucl. Phys. B430 (1994) 577. 

[27] P. Dorey, I. Runkel, R. Tateo and G. Watts, Nucl. Phys. B578 (2000) 85. 

[28] F.C. Alcaraz, M.N. Barber, M.T. Batchelor, R.J. Baxter and G.R.W. Quispel, J. Phys. A20 
(1987) 6397. 

[29] P. de Sa and A. Tsvelik, \cond-matt/ 9 503031 
[30] A.M. Polyakov, J.E.T.P. Lett. 12 (1970) 381. 

[31] A.A. Belavin, A.M. Polyakov and A.B. Zamolodchikov, J. Stat. Phys. 34 (1984) 763; Nucl. 
Phys. B241 (1984) 333. 

[32] H.W.J. Blote, J.L. Cardy, M.P Nightingale, Phys. Rev. Lett. 56 (1986) 742; J.L. Cardy, Nucl. 
Phys. B270 (1986) 186. 

[33] I. Affleck, Phys. Rev. Lett. 56 (1986) 746. 

[34] H.J. de Vega and F. Woyanorovich, J. Phys. A25 (1992) 4499. 

[35] A. Doikou, J. Phys. A36 (2003) 329. 

[36] R.J. Baxter, Ann. Phys. 70 (1972) 193; J. Stat. Phys. 8 (1973) 25; Exactly Solved Models in 
Statistical Mechanics (Academic Press, 1982) 

[37] V.E. Korepin, Theor. Math. Phys. 76 (1980) 165; V.E. Korepin, G. Izergin and N.M. Bogoli- 
ubov, Quantum Inverse Scattering Method, Correlation Functions and Algebraic Bethe Ansatz 
(Cambridge University Press, 1993). 

[38] A. Kirillov and N.Yu Reshetikhin, J. Sov. Math 35 (1986) 2621; A. Kirillov and N.Yu 
Reshetikhin, J. Phys. A20 (1987) 1565. 

[39] L.D. Faddeev and L.A. Takhtajan, Russ. Math. Surv. 34, 11 (1979); L.D. Faddeev and L.A. 
Takhtajan, J. Sov. Math. 24 (1984) 241. 

[40] L.A. Takhtajan, Phys. Lett. A87 (1982) 479. 

[41] A. Doikou and A. Babichenko, Phys. Lett B515 (2001) 220; A. Doikou, Nucl. Phys. B634 
(2002) 591. 

[42] A. Bytsko and A. Doikou, in preparation. 

[43] V.V. Bazhanov and N.Yu. Reshetikhin, Int. J. Mod. Phys. A4 (1989) 115. 



20 



[44] S.R. Aladim and M.J. Martins, J. Phys. A26 (1993) 7287. 



[45] P. Fendley, cond-mat/ 9304031 



[46] D. Bernard, Phys. Lett. B279 (1992) 78. 
[47] N.Yu. Reshetikhin. J. Phys. A24 (1991) 3299. 

[48] G.E. Andrews, R.J. Baxter and P.J. Forrester, J. Stat. Phys. 35 (1984) 193. 
[49] N.Yu Reshetikhin and H. Saleur, Nucl. Phys. B419 (1994) 507. 
[50] E. Date, M. Jimbo, T. Miwa and M. Okado, Lett. Math. Phys. 12 (1986) 209. 
[51] R.I. Nepomechie, \hep-th/02lT00l\ 

[52] A. Doikou and R.I. Nepomechie, J. Phys. A32 (1999) 3663. 

[53] C. Ahn and W.M. Koo, \hep-th/9708080\ J. Phys. A29 (1996) 5845. 

[54] R.E. Behrend, PA. Pearce and D.L. O'Brien, J. Stat. Phys. 84 (1996) 1. 

[55] M.T. Batchelor, V. Fridkin, A. Kuniba and Y.K. Zhou, Phys. Lett B735 (1996) 266. 

[56] P. Baseilhac and K. Koizumi, Nucl. Phys. B649 (2003) 491. 

[57] P. Baseilhac, A. Doikou and K. Koizumi, in preparation. 



21 



